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We study the ryNN couphng constant using the method of QCD sum rules starting from the 
vacuum-to-eta correlation function of the interpolating fields of two nucleons. The matrix element of 
this correlation has been taken with respect to nucleon spinors to avoid unwanted pole contribution. 
The SU(3)-flavor symmetry breaking effects have been accounted for via the ?7-mass, s-quark mass 
and eta decay constant to leading order. Out of the four sum rules obtained by taking the ratios of 
the two sum rules in conjunction with the two sum rules in nucleon mass, three are found to give 
mutually consistent results. We find the SU(3) breaking effects significant, as large as 50% of the 
SU(3) symmetric part. 
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I. INTRODUCTION 



o 

The knowledge of eta-nucleon coupling constant, gr/NN, has implications in both hadronic physics as well as nuclear 
physics. It is required for study of ry production off a nucleon target, and also for analysis of NN scattering data. In 
Qh^ general, it has been used for construction of realistic nuclear potentials [ll-Q- In particular, r?-exchange along with 
^ ^ pion exchange between nucleons can give rise to isospin violation in nuclear potential Q; this has been used for the 
^ resolution of the Nolen-Schiffer anomaly QtiNN will also be useful for understanding the strange content of the 
, nucleon [6|] via the Sullivan process Q • It may provide important constraints on parton distribution amplitudes (DAs) 
of rj. For mesons tt'^, 77 and 77', the strong interaction induces transitions between quarks of different flavors (uu, dd, 
ss) or gluons (gg, ...) The mixing phenomenon is strongly connected with the U(l)yii anomaly of QCD. Hence, it 
is expected that a reliable determination of QtiNN would shed considerable light on the U(l)yi dynamics of QCD 
[■^ — The measurement of gr/NN is a formidable task since the production of rj mesons from single nucleons is dominated by 
. the resonance S'ii(1535) irrespective of the probe Therefore, a reliable theoretical estimate of grjNN is desirable. 
■ Among the various methods used for calculating hadronic parameters, QCD sum rules are especially useful. This 
" I method has already been used for calculating several meson-baryon couplings [lll - l20| and, in particular, grjNN in 
SU(3)-flavor symmetry limit [l5| . In this work we wish to apply this method to calculate grjNN with leading 
, SU(3)-flavor violating effect taken into account. 

' QCD sum rules for calculation of meson-baryon coupling constant was first used by Reinders et al. [Tl], [l^] who first 
" , I ' considered three-point correlation function for g-nNN- On finding of the inconsistency with the Goldberger-Treiman 
J> . (GT) relation they studied two-point correlation function of interpolating fields of nucleons with the pion in the initial 
state. In the former approach one had to assume an extrapolation from large space-like momentum to zero momentum 
^ j for the meson whereas no such assumption was needed in the latter. It was pointed out by Birse and Krippa fl3| that 
in the soft pion limit this latter sum rule can be obtained from nucleon mass sum rule by chiral rotation, and hence it 
does not constitute an independent sum rule from that for the nucleon mass. Therefore, one must take finite meson 
momentum to arrive at an independent sum rule. 

The two-point correlation function with an initial meson state gives rise to several independent sum rules which 
have been studied in detail It has been observed [3 that the sum rule with the tensor structure gives 



the most reliable result. An improvement of this approach has been proposed by Kondo and Morimatsu [18[ wherein 
one takes matrix element of the correlation function with respect to nucleon spinors. The sum rule obtained in this 
way is free from ambiguity arising from the choice of the effective interaction Lagrangian and the effect of nonzero 
mass of the meson can also be taken into account. 

The evaluation of two-point correlation function of nucleons between vacuum and a one meson state requires 
knowledge of matrix elements of nonlocal quark and gluon operators between vacuum and one meson state. Such 
matrix elements have been extensively studied in literature mainly for isospin-non-singlet members of the octet 
pseudo-scalar family [2l| - |23| . There have also been derivations, using various approaches, of Gell-Mann-Okubo type 
of relations involving two-parton and three-parton light-cone distribution functions of pions, kaons and eta .24, .2S]. 
These relations can be used to get matrix elements of nonlocal quark and gluon operators between vacuum and one 
eta-meson state to leading order in SU(3)-breaking corrections. It may be pointed out that constants appearing in 
DAs, such as meson decay constants and others to be specified below, enter both in the effective Lagrangian relevant 
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for low-energy physics and the hght-cone wave functions utihzed in high energy reactions. 

In Sec. II, we derive the sum rule giving necessary details for the projected correlation function. In Sec. Ill, we 
analyze the results numerically and discuss them. Finally, in Sec. IV we summarize our work and give conclusion. 

II. FORMALISM AND CONSTRUCTION OF SUM RULE 

We consider the correlator of the standard nucleon currents between vacuum and one rj-state: 

U{q,p)^t I d^xe'"^ <0\T{Mx),JNm\v{p) >, (1) 
where Jn is the standard proton current 

Jjv = e'"'^[w'^^C7p^i^]757^d^ (2) 
where a, b, c are color indices. The most general form of n(g,p) is [l^ 

u{q,p) = Z75pn^^ + i75n^^ + 75'T^'gpP.n^ + tisgti^^. (3) 

The T^NN coupling constant g^jNN is defined through the coefficient of the pole as [l^ 

u{qr) {q - m) Il{q,p) {q - p - m)u{ks)\q2^^2^ (g-p)2=m2 = iX^griNNu{qr)-f5u{ks), (4) 

where k = q — p and u{qr) is a Dirac spinor with momentum q and spin r and is normalized as 

u{qr)u{qr) — 2m. (5) 

Following [isj. we define the projected correlation function 

n+((7,p) ^ u{q)-f()U{q,p)'yau{q - p). (6) 

n+ can be regarded as a function of qo in the reference frame in which q = 0. We split the projected correlation 
function into even and odd parts as 

Uf{ql) = ^[U{qo)+n{-qo)], (7a) 

n%l)^^mqo)-Ui-qo)]. (7b) 
The dispersion relation satisfied by these functions are 

n^(go') - -- / dqo'^^^ImU+iqo'), (8a) 



We define Borel transform by B 



(8b) 



(-9o) i^ji 



Bf[ql)= fim ^ ' ' , (tt) /(gg). (9) 



so that the dispersion relations become 



ml{ql)] = iy dqo'<zo'ea:p(-^)/™n+(g;), (10a) 
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m%l)] = dq,exp[-^)Imn+{q,). (10b) 

The r.h.s. of Eqs. (|10al) and (jlOb[) are expressed in terms of the observed spectral function. The absorptive part of 
the projected correlation function is parameterized as 



ImU+{q,p) = u{q)i-f5u{q')TrX g{qa,p )[ 



qo - niN 



qo — Ek — UJp 

S{qo - E, -u,)^^ j^^^^ _ ^^^^ ^ ^^^^j ImnOPE{q,p), (11) 



where is the effective continuum threshold of rjN or rjN channel and m^v is the mass of the proton and A is the 
coupling of the proton current with one-proton state: 

< 0| J7v(0)|g Au(g). (12) 

We use the followinR parameterization of the vacuum-to-eta matrix elements of the light-cone operators arising in our 
calculation [3 [H 0, [H {q = u,d): 

< 0|g(0)z75g(x)|77(p) ^[1 - Ip . a; - i(p • x)% (13a) 

/, 1 i 5i 5 1 

< 0|g(0)7p75g(a;)|?7(p) >= -^[^P^^ + ' ^''^^ ^ Ts'^^^ ' 36^^^^^'' ~^ 72^^^ ' ^^'^^'^ ^ 36*^^^^ ' ^''^^^1' ^^^^^ 

< G\gsGl,{]^x)qa{x)qmHp) >= ^C75/(aApP. " fT^.p^) 

+ ^^^^ Cb[g75(7pPt. - luP,,)px + j5P{Pt.x^ - P,^Xf,)]j^-j=f^tef,^°'^[i-/aPp{-^ - -^px)]. (13d) 

The sign of is as per suggestion of Doi et al. [l9j . To leading order in SU(3)-flavor symmetry breakin g, w e use the 
following Gell-Mann-Okubo-like relations for the parameters of octet pseudo-scalar mesons tt, K and rj |23.[25j: 

3/^ + = (14a) 
5fn^^n + I-^y-n = '^Jk^'Kj (14b) 

3/,/i,(l - pT) + - P-) = (1 - P?)- (14c) 

The constants for pions and kaons are taken from Ref. ^]. As a result, we get the following values for the constants 
for ?7-meson: 

= 0.1695Gey, = 1.6380GeV, /s^ = 0.0045Gey^ p[ = 0.1028. (15) 

The decay constant irj obtained is close to the value obtained in a more elaborate scheme of two couplings and two 
mixing angles [26]. The constant /i^ is also close to the value directly obtained for 77 as /i^ = 3TO^/(2m + 4ms) [25| . 
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The expression for the correlatfon function obtained through operator product expansion is: 



+750-'^ qi,p^-^[—-^ln{-q - j ~ 4 +^ 4 < • Gq > 

V6 127r^ 3 firj q^ 27 2^^ 

81 V6 TT q^ 18 v6 TT g^ v6 9 g^ 36 

< —G^ > - Vs- < ^.g^cr -Gqy — + — < —G^ > -F^] + «75^p • q 



' 36g4 ^ TT ^ 18 ^-"^ ^ q^ ■ 81 TT q^' ' '°V6 

r-r 1 7 / 2n fir, _ „ 1 1 ^2 /I 32(^2 _ /i„ _ 1 

96 108 TT g6 9 ^1 



<-qgsa-Gq>-^+J^<^G->- + -<qq>-}]. (16) 



In above expression, q in < qq > and < ggsC • G(7 > stands for u or (i. The above expression has some differences 
from the corresponding expressions obtained previously by other authors [H, [lB| . We have given complete expression 
for the mixed condensate term, and sign of S"^, as pointed out earher, is consistent with Ref. [l^. From Eq. (fT6|) . we 
obtain U^{qQ) and U^{qQ), as given in Eqs. ([5]), ([7a)) and (|7b|) . Using energy-momentum conservation, it can be 
shown that (k ~ q — p) 

2 

771 

f = -77il + {7n^^)/{477i%), c^p = (m^)/(2TOjv), (1^) 
Upon Borel transform, 11^ gives the fohowing sum rule: 

Am.) r 1 ..6 p / s.) ^/-Sfe + - 1 5(5^ (E'fc + 771n ~ ujp) 3 /s,, 21 

" f^^^ ^^^i^^^ — ^^i;^ + ^^'^j^^^'-2 3^1:^ ii^;""^ 

+ {Ek + 771N - i^p){l <qq> -i— < — >)M^ - TTT < — G^ > + {Ek + niN - ujp){^ < qgsCr ■ Gq > 

y y /i^ TT Z4 TT Ot> 

2 Q^s ^9 M /rjA^T) r 9 r / n 7(52M2 4 1 1 Q^s ^9 

< — G^ > + 2^0 + Q < > +T^— < — G' > 

81 TT V6 SflrfTT'^ M 97r-^^,, 9 18 /l^ TT 

+ 18^ < ^9sO- ■ Gq > -^^j;- < ^G^ >} + (Ek + mjv - ojp)pl{^ < qq > + ^ < Q9s<J ■ Gq >} 

+'™'{lT^^o(^) + 4 < — G' >} + ^piEk + mN){-^ <qq> —M^ + ^— < qq > -— < qg^a ■ Gq > 
' 127r'^ 72 TT 3 /Lt^ 27 /i,, /z,, 

-T58 < 

On the l.h.s. / is a constant independent of M"^. On the r.h.s., Eo{x) = 1 — and Ei{x) = 1 — (1 + x)e~^ are 
used to model contributions of excited states. The expression in the second square bracket on the r.h.s. in Eq. (jl8p 
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gives purely SU(3) symmetry breaking contribution. After Borel transform H*^ gives the following sum rule 

2mNX^e~^{gnNN ~ hM^) 

fnH,, J s,, 4 < gg > „ , r ^6 < > ^ 1 < gg > 

12^ + - c..) + ^ - ,,6 ' } < -^G^ > +PliE. + < > Hipl - 1)^ 

+ + - ^p)}^ < ^"'^ > +|[^ < V^' ^ ^""^ " + - (19) 

In Eq. (jl9p also, the expression in the second square bracket on the r.h.s. gives purely SU(3) symmetry breaking 
contribution. On the l.h.s., ft, is a constant independent of Af^. We also write the nucleon mass sum rules [26,,^] for 
our use. The chiral-even mass sum rule is 

m^A^e-^ = ±^[^M^E^{^) < qq > +^ < 99 >< ^G^ >]■ (20) 

In above equation, sq is the continuum threshold for the mass sum rule. The coefficient of dimension 7 operator has 
been taken from Ref. [l^]. The chiral-odd mass sum rule is 

X'e ^ = ^2(^) + —Eoijf,)-' < -G' > +— < m — < ^9s<r • Gg X gg > — ]. (21) 

In Eq. ()2ip . E2{x) = 1 — (1 + a; + x'^ /2)e~^. By dividing the two coupling sum rules with the two mass sum rules 
we get a total of four equations, each of which can be used independently for the determination of the coupling g^NN- 
We will study full sum rules as well as SU(3) symmetric part of sum rules in the next section. 



III. ANALYSIS OF SUM RULES AND RESULTS 

We have used the following numerical constants in our analysis 

< gg -(0.23 Gey)3, < — >= (0.33 Gey)^ < g^.cr • Gg >= < gg > 

TT 

with ml = 0.8Gey^ /a,, = 0.0045, ttiat = 0.939 GeV, m,, = 0.547 GeV, = 0.2 GeV^. (22) 

For continuum threshold s,j in the coupling sum rules we have uniformly used 3^ — 2.57 GeV^ in both coupling 
constant sum rules while so — 2.00 GeV^ in both mass sum rules. Later on we shall study the effect of changes in 
the continuum thresholds on results obtained for g^NN- We use a notation in which f^^^i is the function obtained by 
taking the ratio of full OPE side of even sum rule in coupling to odd sum rule in mass while ^ similar meaning 

except that in coupling sum rule only SU(3) symmetric terms are taken, and similarly for others. We have shown the 
plots of all the functions obtained by taking the ratio of a coupling sum rule to a mass sum rule. Furthermore, for 
coupling sum rules both the full expressions as given in Eqs. and (IT5)) as well as SU(3) symmetric parts have 
been considered. We have found straight line fits for all the curves over a range of Borel mass 1.0 GeV^ < < 1.8 
GeV^. To get a quantitative idea of the goodness of fit, we define a by 

n 

X' = - ff^ti^^)r/{.f{^^) + (^0 }'] / (1 + n). (23) 

i=0 

All the curves obtained from ratios and their straight line fits are shown in Figs. [1] to [121 The Borel windows for both 
coupling sum rules are common, since, in principle, couplings for all the Lorentz structures are related under SU(3) 
rotations [lB| . We have listed the results of fits along with 's in Table HI It is observed that ratios containing even 
mass sum rule, namely, f^f^i and /j^^;, exhibit poor fits. Hence, we have also shown their fits over a shorter Borel 
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TABLE I: Numerical values obtained for g^NN from various sum rules and the (defined in Eq. (|23|) ) for straight line fit over 

the Borel window given in second row. For all the cases n = 20 has been used to get x^- 

Sum rule Even-odd Odd-odd Even-even Odd-even Even-even Odd-even 

Range of (GeV^) 1.0-1.8 1.0-1.8 1.0-1.8 1.0-1.8 1.2-1.8 1.2-1.8 

QnNN in SU(3) sym. limit 3.73 5.31 5.19 9.98 4.83 9.36 

p^jvivfor full expression 5.82 4.20 8.30 6.29 7.60 5.82 

for St. line fit of /toi„, 2.61x 10^^ 1.98x 10"^ 10. 6x 10"=^ 9.23x 10"^ 4.55x 10"^ 4.00x 10"^ 



window ( 1.2 GeV^ < AP < 1.8 GeV^) in Figs. [S] and HUl This gives better x^ but still not so good as that for /°°; 
and f^^cii- Hence, hereafter we analyze these two latter sum rules. 

It is to be noted that contains combination of SU(3) symmetric axial vector and pseudo-scalar terms whereas 
contains tensor terms in SU(3) symmetry limit. The use of tensor structure has been advocated on account of its 
nice features. It has been concluded that for the later the physical parameter, in SU(3) symmetry limit, is independent 
of the form of the nucleon current for currents with no derivative [l5| . the sum rule does not depend on the form 
of the coupling in the effective Lagrangian Borel curves are rather insensitive to the continuum threshold [l5| . 
etc. In the present case, we see from Table |T] that the odd-odd sum rule has the lowest x^- We have found that the 
contribution of excited states to Il^j? is 6.5% at = 1.8% GeV^ while the operators of highest dimension contribute 
6.7% at — 1.0 GeV^. Hence, the Borel window is quite appropriate for this case. The change of continuum 
threshold Sj, from 2.57 GeV^ to 2.07 GeV^ (see Fig. [TT|) increases the value of grjNN by 2.3% while the change of the 
continuum threshold sq from 2.0 GeV^ to 2.3 GeV^ increases it by 3.7%. Change oi < qq > in both the coupling and 
mass sum rules of f^^^i —(0.23 ± 0.025 GeV)'^ changes g^jNN by 9.6%. Uncertainties of gluon condensate, ml and 
6^ may change the value of g^NN by up to 5%. The uncertainties in /x^, /a^ and (which are small); higher 
powers of rriri not considered and extrapolation of the coupling constant to the physical point may result in another 
couple of percent of change in grjNN- In all, we expect an uncertainty of up to 25% in the value of grjNN obtained 
from analysis of /^'^J^; due to uncertainties in various constants used in the sum rule and extrapolation. 

The even-odd sum rule has marginally higher than the odd-odd case. Also, the value of grjNN obtained from 
even-odd sum rule coincides with the one obtained from odd-even sum rule, albeit with a lower Borel window (1.2 — 1.8 
GeV^). Moreover, this sum rule has a minimal dependence on the continuum threshold s,j, as is evident from Fig. 1121 
However, contribution of excited states as per the model used here is 23% at = 1.8 GeV^. Highest dimensional 
operators contribute 24% to H;^ in the entire range of Borel window. This means that the omitted higher order terms 
in OPE may contribute significantly to H;^. Changing < qq> to —(0.255 GeV)'^ in both the coupling and the mass 
sum rule of /j^^; changes g,jNN by 23%. This means that the result obtained from even-odd sum rule has a larger 
uncertainty. Hence, we consider the result obtained from odd-odd sum rule as the most reliable and as our main 
result we have: 



5,,7VAr = 4.20 ± 1.05 (24) 

A few comments on the size of SU(3)-breaking effects in the four sum rules. For even-odd and even-even sum rules, 
the coupling grjNN increases on introducing symmetry breaking effects and the increments (with respect to the SU(3) 
symmetric results) are almost uniform (by 56 to 60%) in these cases. On the other hand, the coupling decreases for 
odd-odd and odd-even cases on introducing symmetry breaking effects. Furthermore, the change is minimum (about 
21%) for odd-odd sum rule and moderate (about 38%) for odd-even case. 

In Table |TT1 we compare our result with the values of g^jv n obtained by other authors using various methods in the 
recent literature. In Refs. [1, |^, [29*1 Goldberger-Treiman relation has been used to calculate grjNN and grj'NN simulta- 
neously. Authors of Refs. [8,9] use the measured values of singlet and octet axial charges of nucleon. Nasrallah |i29| . 
on the other hand, has used dispersion relation for the nucleon matrix element of the divergence of axial current 
saturating it with rj and rj' states and accounting the continuum contribution through a model. Feldmann 's'l has 
neglected the contribution of higher excited pseudo-scalar states. In the absence of knowledge of gluon topological sus- 
ceptibility and gluonic coupling to nucleon (goNN), Shore [9] ends up getting a plot of grjNN and goNN as a function 
of grj'NN ■ Kim et. at. have performed calculation of grjNN to leading order in by considering the two separate 
Dirac structures i^^ljiP^ and and "/scr^i^q^p'^ separately. They have introduced SU(3)-breaking only through the 
decay constant /,,. In Ref. while fitting the data in a potential model that accounts for various meson exchanges, 
the coupling grjNN was not searched independently; rather it was constrained via SU(3) symmetry. In Ref. [3l| the 
authors have worked with effective Lagrangian approach to calculate i] photoproduction in the N*(1535) region. They 
find that the observed differential cross section is not very sensitive to either the nature of the 77-nucleon coupling 
or to the precise value of the coupling constant. They extract a broad range of values for the r]NN pseudo-scalar 
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coupling constant: 0.2 < grjNN < 6.2 from analysis of all available data. 



TABLE II: Some of the results for eta-nucleon coupling obtained in recent literature. 



Ref. no. 


[29] 


f8] 


[14] 


[3] 


[30] 


M 


[31] 


This work 


g-riNN 


5.20±0.25 


3.4±0.5 


5.76 


7.95 


6.8 


3.96 ± 0.16 


3.2±3.0 


4.20± 1.05 








[SU(3) sym ] 






{Qv'nn = 2.0) 












7.34 






3.59±0.15 












[Brok. SU(3)] 






{9,rNN = 1-0) 






Specification 


GT relation 


GT relation 


Sum rule 


Pot. model 


Pot. model 


GT relation 


Expt. 


Sum rule 



IV. SUMMARY AND CONCLUSION 

Eta-nucleon coupling constant has a definite role in analyzing experiments involving NN scattering in general and 77 
production in particular. However, due to the dominant role of N*(I535) in rj production and lack of precise knowledge 
of QriNN, the role of ?7NN vertex is usually ignored in analysis pj, .32] . Hence a credible determination of gr^NN is 
desirable. Separating the even and odd parts in the projected correlation function, we obtained two independent sum 
rules for the coupling griNN- These sum rules incorporate the effect of SU(3)-flavor symmetry breaking in leading 
order. The symmetry breaking parameters themselves come from matrix elements of quark and gluon operators 
between vacuum and one-eta state. Taking the ratios of the two coupling sum rules with the two mass sum rules for 
the nucleon, we obtained a total of four sum rules. We have analyzed the full expression for the four cases as well as 
their SU(3) symmetric forms. The results obtained from three of these full sum rules are quite close to each other and 
this enhances the credibility of our result. The result from odd-odd sum rule is quite robust and is the main result of 
this work. 

We have ignored the mixing of octet state with singlet and glueball states and the consequent role of instantons 
in our analysis. Nevertheless, we have quantified the SU(3) breaking corrections to eta-nucleon coupling constant to 
leading order in SU(3) symmetry breaking, and found the effects to be significant. We conclude that inclusion of 
SU(3) symmetry breaking effects in a realistic way in this work is a step forward in realizing a correct estimate of this 
poorly determined quantity. 
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FIG. 1: Plot of leading order terms in the odd-odd sum rule as a function of (solid curve). A straight line fit of the form 
5.20M2 + 5.31 (broken) over the range 1.0 GeV^ < < 1.8 GeV^ gives gr,NN = 5.31. 




FIG. 2: Plot of full expression obtained in the odd-odd sum rule as a function of (solid curve). A straight line fit of the 
form 1.96M^ -I- 4.20 (broken) over the range 1.0 GeV^ < < 1.8 GeV^ gives g^NN = 4.20. (defined in the text) for this 
fit is 2.0 X 10"^ 
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FIG. 4: Plot of full expression obtained in the odd-even sum rule as a function of (solid curve). A straight line fit of the 
form -1.58M^ + 6.29 (broken curve) over the range 1.0 GeV^ < < 1.8 GeV^ gives g^NN = 6.29. (defined in the text) 
for this fit is 9.2 x 10"^. A better fit is shown in Fig. [H 
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FIG. 8: Plot of leading order terms in the even-evn sum rule as a function of (solid curve). A straight line fit of the form 
-1.42M^ + 5.19 (broken curve) over the range 1.0 GeV^ < < 1.8 GeV^ gives q^nn = 5.19. 
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